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Abstract 

We develop a simplified method for obtaining higher orders in the perturbative expansion 
of the singular term A (q^) /[I — x]^ of non-singlet partonic splitting functions. Our method 
is based on the calculation of eikonal diagrams. The key point is the observation that the 
corresponding cross sections exponentiate in the case of two eikonal lines ||, , and that 
the exponent is directly related to the functions A (a<j) due to the factorization properties of 
parton distribution functions. As examples, we rederive the one- and two-loop coefficients 
and A^^\ We go on to derive the known general formula for the contribution to A^"'^ 
proportional to Nj~ , where Nf denotes the number of flavors. Finally, we determine 
the previously uncalculated term proportional to Nf of the three-loop coefficient A^^^ to 
illustrate the method. Our answer agrees with the existing numerical estimate Q. The exact 
knowledge of the coefficients A^"'^ is important for the resummations of large logarithmic 
corrections due to soft radiation. Although only the singular part of the splitting functions 
is calculable within our method, higher-order computations are much less complex than 
within conventional methods, and even the calculation of A^^^ may be possible. 
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1 Introduction 



1.1 Parton Distribution Functions 

Parton distribution functions (PDFs) are indispensable ingredients in any calculation within 
perturbative quantum chromodynamics involving hadrons in the initial state. fa/A{x) describes 
the distribution of parton a in hadron A with a momentum fraction x. It cannot be calculated 
within the framework of perturbation theory. However, its evolution is perturbatively calculable 
from the renormalization group equation [R| 



(1) 



where Pab is the evolution kernel or splitting function, and fj, denotes the factorization scale, 
usually taken equal to the renormalization scale. This follows from factorization which 
enables us to write a large class of physical cross sections as convolutions of these PDFs with 
perturbatively calculable short-distance functions. Similar considerations apply to partonic cross 
sections. There the parton-in-parton distribution functions //,// describe the probability of 
finding parton /j in parton fj. The evolutionary behavior of the partonic PDFs obeys the 
same equation, (|I|), as for the hadronic PDFs. Thus the splitting functions can be computed in 
perturbation theory. 

The convolutions of these PDFs with short- distance functions simplify to simple products 
when we take moments: 

(iV,^) = = /^"rfxe-^(i-^') (a:, ^) + (^) • (2) 

The first definition in Eq. (|]) is the Mellin transform which is equivalent to the second definition, 
the Laplace transform, at large because then e"^"^^^^') ~ . Tildes here and below denote 
quantities in moment space. 

Up to corrections of order we can neglect flavor mixing, that is, we deal with non-singlet 
distributions: 

The solution to the evolution equation ([l|) for non-singlet parton-in-parton distributions is in 
moment space given by 



ff{N,fx,e) = exp 



12 



/i'2 



■7//(iV,a,(/i')) 



+ (4) 



where 'yff{N) = Jq dxx^~^Pff{x) are the moments of the splitting functions, and where we have 
imposed the boundary condition ff{N,fi = 0,£:) = 1. The argument e in ff indicates that we 
compute in the usual modified minimal subtraction scheme (MS) in n = 4 — 2e dimensions. 
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To leading power in the moments of the partonic sphtting functions, '~fff{N), take the 
simple form 0, 

^ff{N,as) = Afias)\nN + Bf{as) + O (l^) , (5) 



N 



or m x-space, 



1 — X 



+ Bf{a,)6{l - x) + O {[1 - xf) 



(6) 



where the plus distribution 



is defined by 



/ dxf{x) 

J z 



X 



/(I) 



X 



+ /(l)ln(l-^) 



(7) 



The term with the plus distribution represents the cancellation of a single overall infrared diver- 
gence. The coefficient of In N in Eq. (|^) can be expanded in the strong coupling, 



Af{as 



E 



A 



(n) 



The exact knowledge of the terms A^^^ is important for a; — > 1 (large N), since there large 
logarithmic corrections arise due to soft-gluon radiation. These corrections need to be resummed 
in order to be able to make reliable predictions within perturbation theory. The knowledge of the 
coefficients A^^^ and B^'^^ is required for threshold resummation |^ at the next-to-next-to-leading 
logarithmic (NNLL) level 0. 



The anomalous dimensions 7//(A^) are currently known to two loops ||rT|, and a general for- 
^-terms of 7//(-/V) was computed by Gracey [0. Work on the 3-loop splitting 

From the known 



mula for the a^Nj 



functions, based on the operator product expansion (OPE), is in progress [jl3 



exact values for some specific moments and the behavior at small x [|I4[ a numerical parametriza- 
tion for the coefficient A^^^ was obtained by Vogt [Q, although, for the above reasons, the exact 
knowledge of this term is desirable. The above mentioned calculation by Moch, Vermaseren, and 
Vogt of the full three-loop splitting functions via the OPE in moment space will be completed in 
the near future |T^. Their results for the fermionic contributions are now available [|T6|. How- 
ever, the method presented here, although only applicable for the calculation of the coefficients 
A, not of the complete x-dependence, is complementary to the OPE method in two ways: we 
calculate only virtual diagrams, and furthermore, it is much less computationally intensive, thus 
a computation of the four- or even higher loop coefficients may be feasible. 

In the following it is convenient to use light-cone coordinates, where our convention is as 
follows: 



k- 

ki 



-L(.0 + .3), 

{k\e). 



(9) 
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Factorization allows us to define PDFs in terms of nonlocal operators. At leading order in 
one can define P, the following function 



1 

1 



dy e 



dy e 



V 



P 



J2(p\^^{0,y~,0±) n) 7+ 1^,(0)1 p) 



(10) 



which describes the distribution of a quark qi, created by the operator ipi, in a. quark q with 
momentum p. The operators are separated by a light-like distance, and are joined with a path- 
ordered exponential, denoted by P, to achieve gauge-invariance. This exponential describes 
the emission of arbitrarily many gluons of polarization in the plus-direction. A^"^^ is the vector 
potential in the fundamental representation. In the second line we have inserted a complete set 
of states, have used the identity 



and have defined 



(11) 



(12) 



with the light-like vector jS'^ chosen in the minus-direction, = l3± = 0. 

The Feynman rules for the expansion of an ordered exponential are shown in Fig. |l], where 
the double lines represent "eikonal" propagators. Thus, path-ordered exponentials are closely 
connected to the eikonal, or soft approximation, which results in the same Feynman rules, as we 
will see below. Gluon parton distribution functions can be constructed analogously P, |T^, with 
the vector potentials in the adjoint representation, and appropriate operators for the creation of 
gluons. 



1.2 Outline of the Paper 

Below we will outline a method which greatly simplifies the calculations of the collinear coeffi- 
cients in the non-singlet case. Our method is based on the factorized form of the partonic 
non-singlet PDFs at leading order, where for x — > 1 the collinear coefficients are factored into 
an eikonal cross section. It was observed by Sterman that these cross sections exponentiate []1|, 
not only in abelian theories, but also in nonabelian theories. As we will see, this exponentiation 
further simplifies the technical calculations, and enables us to go a step further towards the full 
calculation of the coefficient A^^^ . To illustrate the method, we compute the as yet undetermined 
contribution proportional Nf to the three- loop coefficient. 

A similar observation was made by Korchemsky 0, who related the anomalous dimension of 
PDFs, Eq. (D), with the cusp anomalous dimension of a Wilson loop. His work was performed in 
a noncovariant axial gauge, whereas here we will use Feynman gauge throughout. Korchemsky's 
observation was used in [18| for the calculation of the two-loop coefficient A^'^\ which was done 
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k 
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Figure 1: Feynman rules for eikonal lines in the fundamental representation with velocities (3'^, 
represented by the double lines. The vertical line represents the cut separating the amplitude 
and its complex conjugate. For an eikonal line in the adjoint representation one has to replace 

^ with 2^,. 



in Feynman gauge. The work of Ref. [|TB| was also based on the renormalization properties and 



exponentiation of Wilson loops (see [T^ and references therein). This approach is related to 
ours. However, the additional observations we make below result in several advantages. The 
number of diagrams contributing at each order is decreased by working with light-like eikonals. 
Furthermore, as we will show below, we can restrict ourselves only to virtual graphs. With 
the help of Ward identities we are also able to justify the eikonal approximation and to show 
explicitly the absence of infrared (IR) sub divergences and the cancellation of ultraviolet (UV) 
sub divergences at the eikonal vertex, leaving only the usual QCD UV divergences. 

The discussion consists of two parts. In the first part. Sections |^ and ^, we lay the theoretical 
foundation of our method, whose technical details are described and explained by means of 
various examples in the second part. Sections ^ and ^ 

We will first derive in Sect. Q a factorized form for PDFs, and introduce an eikonal factor 
that absorbs all collinear singularities. This implies that the A{as) can be derived from the 
eikonal factor. In the next section we will review the proof of the exponentiation f^, ^ of this 
factor, show how to construct the exponent order by order in perturbation theory, and prove the 
absence of infrared/coUinear sub divergences at each order. These properties then allow us to 
develop an algorithm, using light-cone ordered perturbation theory (LCOPT), for the systematic 
calculation of the collinear coefficients A, which is equivalent to the calculation of the virtual 
contribution to the anomalous dimension of an eikonal vertex. We summarize the method in 
Section ^, before going on to rederive as examples the one- and two-loop coefficients A^^^ and 
A^'^\ In Section ^ we derive a formula for the coefficients of A^") proportional to which 
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agrees with the corresponding contribution computed by Gracey [|I2| using an effective theory. 
We end by illustrating the steps necessary for the complete calculation of the 3-loop coefficient 
A^^\ The IR structure of A^^^ is explored for the graphs contributing at ce^Nf, which we calculate 
exactly and compare to Vogt's parametrization. In the Appendix we list the Feynman rules 
for LCOPT. 

We will work in Feynman gauge, which further simplifies the expressions and reduces the 
number of diagrams, since then self-energies of light-like eikonal lines vanish. All our explicit 
calculations are performed in the MS scheme, but our method is applicable in any minimal 
subtraction scheme. 



2 Non-Singlet Parton Distribution Functions and Eikonal 
Cross Sections 

Here we will show that the factorized form of a perturbative non-singlet parton-in-parton distri- 
bution function contains an eikonal cross section which absorbs all collinear and infrared singular 
behavior clS X ^ 1 . 



2.1 Leading Regions 

We start with the definition Eq. ([I0| ) of a perturbative parton-in-parton distribution function. 



which is shown in Fig. |^ a) in cut-diagram notation. The part to the left of the cut (the vertical 
line) represents the amplitude, whereas the complex conjugate amplitude is drawn to the right of 
the cut. We pick the incoming momentum to flow in the plus direction, = {p~^, 0, 0_l). Because 
the minus and transverse momenta in (|10D are integrated over, they can flow freely through the 
eikonal line, whereas no plus momentum flows across the cut in the eikonal line. 



We follow the method developed in |21|] to flnd the regions in momentum space which 



can potentially contribute to the leading infrared singular behavior. These regions are shown in 
Fig. 1^ b) in form of a reduced diagram where all off-shell momenta are contracted. Momentum 
regions which can contain singularities are found with the help of the Landau equations p2 



which identify the singularities in Feynman integrals. However, not all regions which satisfy the 
Landau equations give leading contributions. The leading contributions are found by infrared 
power counting pT| . 

We flnd the leading regions for the PDF which are associated with the following momentum 
conflgurations, scaled relative to the large momentum scale in the problem, xp+ ~ p~^: 

• Soft momenta, which scale as A;'^ ~ Xp~^, X <^ 1, in all components. The corresponding 
lines are denoted by S in the flgure. 

• Soft momenta with components which scale in the strongly ordered form a p^,k^ ~ 
Xp~^, X {p^Y, or k" ~ crp+, k^ ~ Xp^, k'j_ X {p'^)'^, respectively, where a <^ A ^ 1. 
These are the so-called Glauber- or Coulomb momenta UM. 
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Figure 2: a) Graphical representation of a parton-in-parton distribution function, Eq. (|TOD, and 
b) its factorized form for arbitrary momentum fraction x, drawn as a reduced diagram, c) The 
reduced diagram for the PDF in the hmit a; — > 1. Here we have suppressed the labels L and R 
for purely virtual contributions to the left and to the right of the cut, respectively, compared to 
the notation in the text. The cut represents the final state. 
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• Momenta collinear to the momenta of initial or final state particles (including the eikonal). 
The corresponding subdiagrams are denoted as J^, where m labels the external mo- 
mentum. Momenta collinear to particles moving in the plus direction scale as ~ 
p'^,k~ ~ Xp~^,k'j_ ~ A (p"^)^, whereas momenta collinear to the minus direction behave 
as k~ ~ p+, /c+ ~ Ap+, ~ A (p"*")^. 

• These regions are connected by a hard scattering, H, where all momenta are far off-shell, 
and thus scale as ~ p"*" in all components. 

In our case we can have a jet Jp collinear to the incoming momentum p~^, as well as an arbitrary 
number of jets Jj emerging from the hard scattering. Furthermore, we can have momenta collinear 
to the eikonal moving in the minus direction, represented by in the figure. The jets can be 
connected by arbitrarily many soft gluons, S. Here and below, unless explicitly stated otherwise, 
we use the term "soft" for both soft and Glauber momenta. 

A further simplification of the leading behavior occurs when we perform the limit to a; ^ 1, 
in which we are interested here. Jet lines having a finite amount of plus momentum in the final 
state become soft, and only virtual contributions can have large plus momenta. This does not 
affect the jet collinear to the eikonal, since it is moving in the minus direction. Thus we arrive 
at the leading regions depicted in Fig. ^ c), with hard scatterings and iJ/j, which are the 
only vertices where finite amounts of momentum can be transferred, virtual jets Jp^L and Jp^R 
collinear to the incoming momentum p"*", a jet Jg collinear to the eikonal (3~, connected via soft 
momenta, S. Here and below the subscripts L and R, respectively, indicate that the momenta 
and functions are purely virtual, located to the left or to the right of the cut. 



Infrared power counting pO|, ^ shows now that the parton distribution function for x ^ 1, 
as shown in Fig. ^ c), can at most diverge as one inverse power, which corresponds to the 
divergence we are expecting from Eq. (j^). Since the power counting is straightforward we just 
state the result here. To get the maximum degree of IR divergence the following conditions have 
to be fulfilled (see the figure): 

• No soft vectors directly attach Hl or with S. 

• The jets and the soft part can only be connected through soft gluons, denoted by the sets 
{qL,mJ, {qR,ma}, and {qi}. 

• Arbitrarily many scalar-polarized gluons, {kL,ij^}, {kR^ij^}, {kLj^} and {k^jj^} , attach 
the jets with and Hr, respectively. The barred momenta are associated with Jg, the 
unbarred with the Jp's. 

• Exactly one scalar, fermion, or physically polarized gluon with momentum kj^ — J2 k^ , 
kf^ — J^kji i^, k^ — J^kl jj^ or A;^ — J2kji j^, respectively, connects each of the jets with 

iR ' jh ' Jr ' _ _ 

the hard parts. The momenta fc^^ {k'ji ) and fc^^ (^^^) denote the total momenta flowing 
into Hl {Hr) from Jp ^ {Jp,r) and Jg, respectively, that is, they are the sum of the scalar, 
fermion or physically polarized gluon momenta, and the scalar-polarized gluon momenta. 
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In an individual diagram we can have only scalar-polarized gluons connecting the jets with 
the hard parts, and no scalar, fermion or physically polarized gluon. However, the sum of 
these configurations vanishes after application of the Ward identity shown in Fig. ^ a). 

• The number of soft and scalar-polarized vector lines emerging from a particular jet is less 
or equal to the number of 3-point vertices in that jet. 

In summary, we have found that the regions in momentum space which give leading contri- 
butions may be represented as: 

d^kr, r d^kR f d^kL f d^kR ^ fd^kL,i^d"kR^,^ d^kL,,^d^kR^,^ 



^ ^^^i (2vr)«y (27r)"i (27r)"i (27r)^^_,f j^^^ (27r)" (27r)" (27r)" (2vr)« 

-pr f d''qL,mL d"qR,nin d'^qi ^cs) ([.ri^L \.i,P^R X-in^iX" 

X H, {kl\ {k^^J; fcf , {kl)j) Hr (fc-', {fc^- }; J) 



X J,,^ {kl\ {k^^J; {ql:^j) J,,R {fc^- }; {q'£-^]] 

X jf^^ (fcf , {fc^^^.J; {kl-^}- {gf'}) (13) 

X 5- (e m.+H q'r, ™« + E ^i] ^" (^L + ki- xp^) 6- {¥r + ¥r- xf 
X 5- [kl - E qU, - 5- (k-R - E q^m, - P"-] ■ 



The sum in Eq. (|T3|) runs over all cuts of jet Jg, C^, and of the soft function S, Cs, which are 
consistent with the constraints from the delta-functions due to momentum conservation. The 
functions in (|13|) are still connected with each other by scalar polarized or soft gluons. This 
obscures the independent evolution of the functions. In the next subsection we will show how to 
simplify this result. 

2.2 Factorized Form of the Perturbative Distribution Function 

Following IP , we will show here that the scalar polarized gluons decouple via the help of a Ward 
identity, and that we can disentangle the jets and the soft part, which are connected by soft 
gluon exchanges, via the so-called soft approximation. 

2.2.1 Decoupling of the Hard Part 

Starting from Eq. (0), we use the fact that the leading contributions come from regions where 
the gluons carrying momenta {A;L,j^}, {kR^i^}, {^1,]^} and {k-R^jj^} are scalar polarized. Thus 
i^i, Hr, Jp,L, Jp,R, and Jg have the following structure: 
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where, as above, for the functions to the right of the cut, we replace the subscripts L with R. In 
these relations the vectors 

e = St, 

= S'l (17) 

are the light-like vectors parallel to p'^ and parallel to the direction of the eikonal, respectively. 

We now use the identity depicted in Fig. | c), where the grey blob denotes the hard part. 
Fig. ^ c) follows from the Ward identity shown in Fig. ^ a), and the identity for scalar polarized 
gluons attaching to an eikonal line in Fig. ^ b). The Ward identity says that the sum of all 
possible attachments of a scalar-polarized gluon to a matrix element vanishes. From this follows 
Fig. ^a), since, by definition, we do not include the graph into the hard function where the gluon 
attaches to the physically polarized parton (quark or gluon), shown on the right-hand side of Fig. 
^ a). The eikonal identity in Fig. ^ b) follows trivially from the eikonal Feynman rules in Fig. 
^ Thus, since the right-hand sides of figures a) and b) are the same (the "empty" eikonal line 
carries no momentum), the left-hand sides are the same. Repeated application of this identity 
results in Fig. |^ c) ||25[. Note that the color factors are included in the Ward identity, resulting 
in the appropriate color factor for the attachments of the gluons as shown in Fig. c). 

Using the identity in Fig. ^ c) for all scalar polarized gluons in the sets {ki^ij^}, 
{^ljl} {kRjj^}, we achieve our goal of decoupling the jet functions from the hard function. 
This decoupling occurs in Feynman gauge only after summation over the full gauge-invariant set 
of graphs which contribute to the reduced diagram Fig. ^ c). The result is shown in Fig. ^ The 
products over the vectors j3 and ^ are replaced by eikonal factors £, which we can group with 
the jets. Furthermore, the hard scatterings, by definition far off-shell, become independent of x 
up to corrections which vanish for x ^ 1. Eq. (^) then becomes 

' d"-(lL, rriL dP'QR, rriR dPqi 



X 



Jdyjdz {yp,f,; {ql-^J; {g^T^J; {gf }) jf'\zp,^-^- {gf }) (18) 

\mL niR I / 



5{l - X - y - z) , 

where /i denotes the renormalization scale, which we set equal to the factorization scale, for 
simplicity. Corrections are subleading by a power of 1 — x. We define the functions Jp,L, Jp,R, 
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Figure 3: a) Ward identity for a scalar polarized gluon. b) Identity for a single longitudinally 
polarized gluon attaching to an eikonal line, c) Resulting identity after iterative application of 
Figs, a) and b). Repeated gauge-group indices are summed over. 



and Jf3 as follows: 



(19) 



and Jp^R is defined analogously to Jp,L, with the subscripts L replaced by R, and with a complex 
conjugate eikonal line, since it is to the right of the cut. In ([T8| ) the total plus momentum 
flowing across the cut is restricted to be (1 — x)p^ , and flows through the soft function and/or 
the eikonal jet. The plus momenta flowing across the cuts Cs and Cp, denoted by yp'^ and zp'^ , 
respectively, are therefore restricted to be (1 — x)p'^ via the delta-function in Eq. ([TSp. Above 
we have factorized the hard part from the remaining functions, which are still linked via soft 
momenta. 



2.2.2 Fully Factorized Form 



Here we will use the soft approximation [24, p5| 



to factorize the jets Jp from the soft 



function, which, by power counting, are connected only through soft gluons. We could factorize 
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the eikonal jet from S in an analogous way, but we choose not to do so here because eventually 
we will combine all soft and eikonal functions to form an eikonal cross section. 

The soft approximation consists of two approximations: the neglect of non-scalar gluon po- 
larizations, and the neglect of soft momenta q'^ in the numerator compared to jet-momenta /c^ 
and fc^ compared to g ■ in the denominator. After making these approximations we proceed in 
a fashion similar to the previous subsection. 

Let us start by decomposing each soft gluon propagator D^i, in the sets {gL,mi,} and {qR^mn} 
into a scalar-polarized contribution and a remainder [p7|] : 

Df,^{qL,mL) = G^^{qL,mL,0 + Kf^^iQL,mL,0 ^ (21) 

and analogously for the momenta {qR,mR} to the right of the cut, where we define 



The sign of the ie-prescription is chosen in such a way as to not introduce new pinch singularities 
near the ones produced by the soft gluons under consideration. For momenta to the right of the 
cut in initial-state jets we also have +ie, since the sign of the momentum flow to the right of the 
cut is reversed relative to the jet momentum k. That is, we have for momenta to the left of the 
cut (g + + i€ ~ 2q- k + ie, whereas to the right of the cut we obtain {q — k)^ — ie ~ —2q ■ k — ie. 

From Eq. (^2]) we see, that for the scalar-polarized ii"-gluons the identity shown in Fig. |^ 
c) is immediately applicable, leading to the desired factorized form. So it remains to be shown 
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that the G-gluons do not give leading contributions. As mentioned in the previous subsection 
power-counting shows that only 3-point vertices are relevant for the coupling of soft gluons to 
jets. Let us consider a 3-point vertex, where a soft G-gluon with propagator Gfj,u{q) couples to a 
fermion jet-line with momentum fc^ in the jet Jp moving in the plus-direction: 

-YjT^G,, (g, ^ 2^^— ^— Ly^-T^'^^M- ^ 0, (23) 



{k + qY + ie' k'^ + ie '^^'^'^^ ^(^k + q)^ + ie) {k^ + ie) 

because k'^ ^ k'^^^ with corrections proportional to Xp^^X <^ 1, as follows from the power 
counting described in Sect. |2TT|, and because G^y {q,0^'^ = 0- analogous observation holds for 



the coupling of G-gluons to jet-lines via triple-gluon vertices. In (p3D we neglect all terms of order 
Xp~^ in the numerator, including the momentum q, because we assume that the denominator also 
scales as ~ Xp'^. This approximation is only valid for soft gluons not in the Glauber or Coulomb 
region, where the denominator behaves as ~ X'^p~^. The remaining step in order to obtain a fully 
factorized form of the PDFs is therefore to show that the soft momenta are not pinched in this 
Glauber /Coulomb region. Then we can deform the integration contours over these momenta 
away from this dangerous region, into a purely soft region where the above approximations are 
applicable. 

Consider again the 3-point vertex in Eq. (p3D, where now the gluon with momentum q is in 
the Glauber /Coulomb region. If is not dominant over \2k± ■ q± + q±\ in the denominator 

our approximation fails. The poles of the participating denominators are in the q~ complex plane 
at 

ql-ie 

q- = (24) 

As long as the jet-line k carries positive plus momentum, we see that the g~-poles are not pinched 
in the Glauber region. In this case we can deform the contour away from this region into the 
purely soft region, where jA;^?^! ^ \2k±_ ■ q^ + qV\- In a reduced diagram, which represents a 
physical process, there must be at every vertex at least one line which flows into the vertex, 
and at least one line which flows out of the vertex. Thus, at every vertex, we can always find 
a momentum k for which the above observation holds, and we can always choose the flow of 
q through the jet to the hard scattering along such lines. Because the soft gluon momenta 
{^L.mi,}, which we want to decouple, connect only to the purely virtual, initial-state jet JpQ, this 
observation remains true throughout the jet. An analogous argument applies to the right of the 
cut. 

In summary, soft gluons can be decoupled from initial state jets at leading power in 1 — x. 
We therefore arrive at the factorized form of the parton distribution function as x — 1: 

/j^i (x) = HL{p,fi;(3,0 HR{p,^I■,(3,^)Jp,L{p,^^■,^3■,0 Jp,R{p,f^;P;0 



^Initial and final states are defined with respect to the hard scattering. 
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X j dy J dz [yp, {gl;- J; {gj^- J; {gf'}) ^^''^(zp, /x; {gf'}) 

X ^"fEC.+E9K,™« + E?T) 5(1-^-1/-^), (25) 

Vm^ TTifl I J 

where we have grouped the eikonal factors stemming from the soft approximation with the soft 
function and the eikonal jet. We define 

^p,L(p,/^;/?;{gZ7mi}) = ^(^,{?L,mi})^^'"^^ ^ L (p, , (26) 

and analogously for the jet to the right of the cut, Jp,_R, with a complex conjugate eikonal. 
2.2.3 Fully Factorized Form with an Eikonal Cross Section 



,{7m„} 



Although in Eq. ( Pol) the various functions are clearly separated in their momentum dependence, 
the parton distribution function is not quite in the desired form yet. We want to write the PDF 
in terms of a color singlet eikonal cross section, built from ordered exponentials: 



(cik) 



^1 — x)p^ 



dy- 



J{l-x)p+y 



Tr 1 7 271 
X Tr (0 



W^'"'\0,y-,0_ 



W^^^HO)! O), (27) 



where the product of two non-Abelian phase operators (Wilson lines) in the representation a, for 
quarks, is defined as follows: 



x] 



<^(A2,Ai;x) 



Pe 



-ig J^^ dA/3-yt(^ )(A/3+x) 



(28) 
(29) 



where the light-like velocities ^ and P are defined in ([ITD , and where A^^^ is the vector potential 
in the representation of a parton with flavor /. The trace in (|27| ) is over color indices. The 
lowest order of the eikonal cross section is normalized to 5(1 — a;). This eikonal cross section has 
ultraviolet divergences which have to be renormalized, as indicated by the renormalization scale 
II. Furthermore, the delta-function for the soft momenta in Eq. (|25|) constrains the momentum 
of the final state in a^°''^^ to be (1 — x)p^. 



We can factorize the eikonal cross section (p7|) in a manner analogous to the full parton 
distribution function, and obtain 



(oik) 



X] 



_ sr^ Y\ f d^(lL,mi^ d'^qR^ruR d'^qi 
~ ctisJLi^ (2vr)" (2vr)" (2vr)- 
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X ^''(T.€m, + T.<n., + T.<ii] ^i^-^-y-^)- (30) 

The eikonal jets Jp'i\ Jp'ji moving collinear to the momentum p, are defined analogously to Eq. 
([T9|), with the fermion line carrying momentum p replaced by an eikonal line in representation a 
with velocity /?. We can define analogous to Eq. (|26| ) 

4:? = ^(e,te,„j)^""^^ ^?(/i;/5;0, (3i) 

and similarly for the jet to the right of the cut, Jp'^ , with a complex conjugate eikonal. In the 
following we suppress the index a for better readability. 

Combining Eqs. (pS]), (pO|), and (pT]), we arrive at the final form of the factorized parton 
distribution function, shown in Fig. ||, 

fp' (x) = Hl {p, /i) Hn {p, /i) Jp, /i) Jp«^(p, /i) ((1 - x)p, fi) , (32) 

suppressing the dependence on the lightlike vectors (|l3). The purely virtual jet-remainders are 
defined by 

The contributions to the various functions in Eq. (|32|) can be extracted from any Feynman 
diagram by examining each region of integration separately. In the following, we only need the 
eikonal cross section in the form (pTl), as well as the x-dependence and renormalization properties 
of the PDF and the terms constituting its factorized form. These properties are the subject of 
the next subsection. 



2.3 Renormalization of Parton Distribution Functions 



As was shown in the parton distribution functions, defined in their unrenormalized form in 
terms of nonlocal operators (0), obey the evolution equation (|^), where the kernel Pat is found 
from the usual relation |T^, ^ 



X 



X 



(34) 



where g^/ (47r) = as, and where In denotes the --pole of the counterterm which multiplies the 
plus-distribution, plus scheme dependent constants, if we work in a minimal subtraction scheme 
with dimensional regularization. Above we only exhibit the term which is singular as x ^ 1, 
since it is this term which we want to extract from the renormalization of our factorized form, 
Eq. m- 
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Figure 5: Parton distribution function for x — > 1, factorized into hard scatterings, an eikonal 
cross section, and purely virtual jet-remainders, as derived in Eq. (p2D. The virtual jet functions 
are normalized by their eikonal analogs, as in Eq. (PBD- 

From Eq. (|3^ we observe that only the eikonal cross section can contribute to the A-teim 
proportional to a plus-distribution. This is because the hard functions are off-shell by O {xp'^), 
and the jet-remainders are purely virtual, thus cannot contain plus-distributions. Therefore, 
their renormalization has to be proportional to Bf 6{1 — x), as was observed in [^]. 

It is thus the renormalization of a color singlet eikonal vertex which we have to study, in order 
to compute the coefficients A^"'^ in (^). Although there are, aside from the usual QCD divergences, 
additional divergences at the eikonal vertex, this is a significant simplification compared to the 
investigation of the UV behavior of the full PDF. Let us now study the eikonal cross section as 
defined in Eq. (|27D, whose perturbative expansion follows the Feynman rules shown in Fig. |1|. 

3 Nonabelian Eikonal Exponentiation 

Although the eikonal approximation simplifies the perturbative calculation already significantly, 
there are still many diagrams to be calculated at each order. In addition, the eikonal approxi- 
mation also introduces new infrared divergences. 

An observation first made by Sterman then proved by Gatheral 0, and Frenkel and 
Taylor solves both problems. This theorem states that a cross section X with two eikonal 
lines in a nonabelian theory exponentiates, 

^(eik) ^ X = e^, (35) 

where Y can be given a simple recursive definition. In this section we will recall the proof of Eq. 
dH) [B @]' f*^^ sake of completeness, including a few illustrative examples which will be used 
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below for the calculation of the 2- loop coefficient A^'^^ . The exponent Y in (|35| ) has the following 
properties: 

1. y is a subset of the diagrams contributing to X, which we will call "webs" in the following, 
since, as we will see below, their lines " [...] are all nested [...] in a spider's web pattern" 

i- 

2. The color weights of the diagrams in Y are in general different from those in X. 

3. For Eq. (|35D to hold, the phase-space region should be symmetric in the real gluon mo- 
menta. 

Below we will outline the arguments necessary to prove this theorem. The proof relies on 
the recursive definition of color-weights and on the iterated application of a well-known eikonal 
identity. Then we go on to show that IR and UV sub divergences cancel in the exponent at each 
order. 



3.1 Proof of Exponentiation 
3.1.1 Some Terminology 

In order to specify which subset of diagrams of the original perturbation series X contributes to 
the exponent Y we need to introduce some terminology. 

Each diagram will be decomposed into its color part and its Feynman integral in the eikonal 
approximation. The eikonal Feynman rules were given above in Fig. |I|. The color part can be 
represented graphically in a diagram which is similar to an ordinary Feynman diagram, but the 
vertices represent the color part of the Feynman rules, i.e. the vertices are just the T-"s and 
i fijkS, for quark or gluon, respectively, and the lines are 6ijS. In addition, all soft lines have to 
be drawn inside the (cut) eikonal loop for reasons which will become clear shortly. Certain color 
diagrams are related to each other by use of the commutation relations of the T^^s and i fijkS 
(Jacobi identity) which are graphically represented in Fig. ^. 



tfabcT" 



filmfmjk ~l~ fjlmfimk ~l~ fklmfijm 0- ("^6) 

As mentioned above, the diagrams contributing to Y will be called "webs" [Q. Originally 
0, a web was defined as a set of gluon lines which cannot be partitioned without cutting at 
least one of its lines. As already stated above, all soft lines are to be drawn inside the eikonal 
loop(s). However, at O [a^) new types of diagrams arise which Frenkel and Taylor called 
connected webs ("c-webs"). c-webs are not included in the original definition for the following 
reason: If one cuts the horizontal gluon line of the c-web drawn in Fig. ^ one would get two webs 
consisting of three-point vertices since real and virtual gluon lines are treated on equal footing 
in a color-weight diagram. Below we will refer to webs and c-webs just as webs. 
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a b a b a b 




Figure 6: Graphical representations of the commutation relations between a) T^jS, and b) i fij^s 
(Jacobi identity), Eq. (|36D . 



The definitions given by Gatheral, and Frenkel and Taylor can be unified by the following 
definition: 

A web is a (sub) diagram consisting of soft gluon lines connecting two eikonal lines which 
cannot be partitioned into webs of lower order by cutting all eikonal lines exactly once. Stated 
differently, webs are two-eikonal irreducible diagrams. The order of a web is defined to be equal 
to the powers of it contains, e. g. a web of O (a^) will be called a web of order 2. Diagrammatic 
examples are shown in Fig. |^. A web has a color factor C and a Feynman integral part JF. JF 
contains only those eikonal propagators which are internal to the web. The color weight is in 
general different from the one which one would get from the usual Feynman rules. 




Figure 7: Examples of a) a web (order 2), b) a c-web (web of order 3), and c) a diagram which 
is not a web (consisting of a product of an order-1 web and of an order-2 web). 

The color weight of a web of order m is recursively defined as 
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^ T^C{W(% (37) 

where CiW^"^^) is the ordinary color factor, is the usual normalization (see Eq. (|27|)), 
C(°) = Tr 1, Ed is the sum over the set of all non-trivial decompositions d of W^"^^ into webs of 
order i < m, and denotes the product of all webs rii of order i (1 < i < m) in a particular 

rii 

decomposition d. The set of all non-trivial decompositions of a given web can be obtained by 
successively disentangling crossed gluon lines in the web by repeated application of the color 
identities given in Fig. ^ 

In [H Gatheral showed that webs in the original definition have what he called "maximally 
nonabelian" color weights ~ a^CpC^'^ , where the CjS are the Casimir factors in the funda- 
mental and adjoint representation, respectively. This statement, however, is misleading at orders 
> 0. We will see an example below, in the calculation of the Nf term contributing to the 
coefficient A at three loops. 

In the following subsection we will clarify the above definitions in an example which shows 
how to factorize eikonal Feynman diagrams into sums of products of webs. This then leads 
directly to exponentiation. The recursive definition of the color weights of the webs ensures the 
factorization of the color parts. For the factorization of the Feynman eikonal integrals we will 
make repeated use of the eikonal identity 



11111 

(3^ 



p ■ kip ■ {ki + k2) p ■ k2P ■ {ki + k2) p-kip-k^ 

illustrated in Fig. |^ a). This identity can be extended to an arbitrary number of soft gluons 
in a straightforward way by repeated application of Eq. (|38D : For two webs Wi and W2 with 
gluon legs ki{i = 1, . . . , m) and Ij {j = 1, . . . , n) attached to an eikonal line with velocity p the 
generalized identity reads 

1 1 

T(W^)J^(W2) ~ 

p-kip-{ki + k2)...p-{ki + ... + kJp-kp-{k + l2)...p-ili + ... + ln) 
- F, (39) 

perms(n,m) 

where the sum is over all Feynman diagrams F obtained by permuting the n + m gluon lines 
such that the order of the ki, and Ij, respectively, within each web is not changed. A simple 
example is shown in Fig. || b). The extension to more than two webs follows by repeating the 



above argument: 



T.F = 1[^(WS). (40) 



Fin d 



3.1.2 An Example 

We will show by induction that the terms in the perturbation series X in Eq. (0), normalized 
by the zeroth order contribution, can be reorganized into a sum of products of webs which can 
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b) 



ki k2 



ki k'l h ki li k2 h ki k2 



Figure 8: a) Eikonal identity for 2 gluons, and b) an illustration of the generalized eikonal identity 
Eq. (|39D for two webs. 



be rewritten as exp(y). Therefore it is necessary and also instructive to start with the first 
nontrivial example, diagrams of 0{a1). At this order we have the Feynman diagrams shown in 
Fig. ^ for quark or antiquark eikonal lines, excluding eikonal and gluon self-energies. Eikonal 
self-energies vanish if we work in Feynman gauge. The sum of diagrams at a given order in is 
gauge invariant, of course. The contribution of the first two terms in Fig. |^ can be rearranged 
by applying Eqs. (^) and (|38D as shown in Fig. |10|. 

Thus we arrive at the following series obtained by rearranging the expansion of X up to 

all webs 
of order 1 



2! . 

all webs all webs 

of order 1 of order 2 



which is illustrated in Fig. ITT]. The combinatorical factor 4 is necessary to avoid overcount 



ing since two webs with the same structure are indistinguishable if the integration measure is 
symmetric in the real gluon momenta. 

3.1.3 Exponentiation 

Looking at the above example it is now clear that any Feynman diagram with two eikonal lines 
can be expressed as a sum of products of webs by applying Eqs. (^) and (^0|) repeatedly. By 
induction we arrive at the following equation for the set of Feynman diagrams of 0{a^), F^*^^: 



1 / \ Ui 

^^"^ = En^ E CiW^'^)J^iW^'^)) , (42) 



{n^} i all webs 

of order 1 
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Figure 9: Diagrams contributing at (excluding self-energies). The color factors are given 

for eikonal lines in the fundamental representation, omitting the overall normalization ev- 
erywhere. 




[c'(im<^)f + c'(X)^i 





b) 

^(X) = - cilf 
Oil I) = c{l I) - Cilf = 



Figure 10: a) Rearrangement of the first two terms of Fig. |^ using Eq. (pHD, and b) Eq. (|37| ) 
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Figure 11: Graphical expression of the exponentiation up to webs of order 2 in Feynman gauge. 
The color weights are given for (anti) quark eikonal lines. 



where i labels the order of the webs and the sum is over all sets {rii} , < rij < oo such that 

Y^iiUi = n. For example, at 0{a^) wc can have ni = 3 webs of order 1 ({3, 0, ...}), orni = 1 webs 
of order 1 and n2 = 1 webs of order 2 ({1,1,0,...}), or = 1 webs of order 3 ({0, 0, 1,0,.. .}). 
The combinatorical factor of is needed to avoid overcounting because of property 3.) of X, 
stated in the introduction to this section, namely that the integration measure is symmetric in the 
real gluon momenta, for example J^^^\ki)J^^'^\k2,ks) — T^^\k2)T^'^\ki,k^), which means that 
webs of the same structure are indistinguishable. Were property 3.) not fulfilled, the perturbation 
series would not exponentiate. 

We now rearrange the original perturbation series given in powers of a" 



oo 



(43) 



n=0 



^ = EE 



n=0 {m} 



i 'H- ^ all webs ^ 



of order t 



all possible i 



{rii} 



E 





(44) 



of order i 
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where we have used the fact that for any function f{ni,i) 

E n/(^-^)=nE/(^-^) (45) 

all possible i i fT'i 

{rii} 

which is easy to see by comparing the expansions of the left and the right hand sides. 
So the series exponentiates 

X = e^, 1^ = E( E C{W^'^)J^{W^'^)]. (46) 

t all webs 

of order i 

This completes the proof that eikonal cross sections with two eikonal lines can be written as an 
exponent of an infinite sum of webs. 



3.2 Cancellation of Subdivergences in the Exponent 



Gatheral, Frenkel, and Taylor showed in by explicit fixed-order calculations that infrared/ 
coUinear subdivergences cancel in the exponent. Here we will outline the proof of this cancel- 
lation, as well as of the cancellation of UV subdivergences involving the eikonal vertex, to all 
orders with the help of the identities in Fig. ^ in the soft approximation. The remaining UV sub- 
divergences are removed via ordinary QCD counterterms, and thus an additional investigation 
of the renormalizability of the eikonal vertex is unnecessary. 

To show the absence of subdivergences, let us rewrite Eq. (^2|) as 

E 

order n {ni},i<n i ^ all webs 



^W^W = + irW _ ^ n— E C{W^'^)J-{W^^^)) . (47) 



of order 2<n 



Eq. { ^7\ j means, that the sum of all webs at order n are given by the original perturbation series 
at that order where all lower-order webs have been subtracted out. The original perturbation 
series can be classified into terms without subdivergences, denoted by Fj:^]^, and terms which 
contain subdivergences, f]"''. 

In eikonal cross sections, infrared/collinear divergences stem from the same momentum con- 
figuration as UV divergences. Since eikonal cross sections are scaleless, when a line becomes 
coUinear to an eikonal it can carry infinite momentum in a light-like direction. But in this case 



we can employ the soft approximation described in Sect. |2.2.2| to factorize these jet-like config- 



urations from the rest of the eikonal cross section. The reasoning follows Section |2.2.2| , and we 



arrive at the equality shown in Fig. |T2| The grey oval in the figure represents a specific jet-like 
configuration, coUinear to one of the eikonal lines. The displayed equality states that the sum of 
all webs at a given order, where this jet-like configuration is connected to the rest of the eikonal 
cross section by soft gluons, can be expressed in the factored form shown on the right-hand side. 
As in Section |2.2.2| , remainders are non-leading. Due to the definition of the color weights (0), 



the right-hand side does not constitute a web of the same order, but rather a product of webs of 
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Figure 12: Factorization of jet-like collinear configurations, represented by the grey oval, from 
the eikonal cross section with the help of the soft approximation. 



lower orders. The contribution shown on the left-hand side of Fig. |T2| is a contribution to f]"'' of 



Eq. (^71). In (P7D, however, we subtract out all products of webs of lower orders, thus cancelling 



the divergent contributions because of the equality shown in Fig. 12. Using the equality in Fig. 



12] and Eq. (^) recursively for every IR/coUinear sub divergence, we see that the sum of webs at 
a given order is free of such sub divergences. 

To summarize, the collinear configuration does not contribute at the order under consideration 
after summing over all relevant webs at that order, because this collinear configuration has already 
been taken into account at a lower order. The only possible collinear and UV vertex divergence 
can occur in the final, overall integral. Of course, in the original perturbative expansion X of the 
eikonal cross section in Eq. (^) these collinear and UV subdivergent configurations contribute, 
but in the exponent Y of Eq. p6| ) they only appear as overall divergences. 

4 The Method for Obtaining Af^ 
4.1 Renormalization of Webs 

As was observed in the above sections, webs have at most one collinear/IR divergence, and one 
overall UV divergence from integration over the final momentum. Using the invariance of eikonal 
cross sections under rescalings of the eikonal velocities, and the fact that there is at most one 
overall IR divergence, we can deduce the dependence of the webs on its overall momentum k, 
and on the light-like eikonal momenta. For a color singlet eikonal cross section we have 



after integration over internal momenta, in a frame where the eikonal lines are back-to-back P8 |. 
For definiteness, we pick the incoming line ^ moving in the plus direction, and the outgoing 
eikonal (3 in the minus direction, and since quantities built from eikonal lines are scaleless, we 
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can scale the eikonal velocities to 1. Recall (|T^: 



e = (1,0, Ox) 

(3 = (0,1, Ox) (49) 

in light-cone coordinates. This choice will simplify the calculations considerably, as we will see 
below. We can therefore write the contributions from virtual webs of order n to the eikonal cross 
section as 



cr'i^i/ri^fml^]^^ ,50) 



where / is an integer < n, and K contains numerical factors (including factors of tt) and is, in 
general, a function of e due to the regulation of infrared and UV (sub) divergences. Above, on 
the left hand side, all internal momenta have been integrated over, as well as and internal 
UV divergences have been renormalized. The integration over fc^ results in terms ~ -^pymF- 
For graphs including (local) counterterms / < n, whereas for graphs with n loops / = n. Both 
virtual webs and their complex conjugates contribute to the overall factor of 2. The structure 
of the integral over follows from boost invariance. In Eq. (^0]), this integral is divergent, but 
these divergences cancel against the corresponding real contributions, and therefore do not affect 
the anomalous dimension of the eikonal vertex. The fc^-integral plays the role of jz^ for the full 
parton-in-parton distribution functions (cf. Eq. (P?])). 

The final scaleless k± integral provides the n-loop UV counterterm which contributes to the 
anomalous dimension Pjj, Eq. (|^). To isolate the UV pole we temporarily introduce a mass 



(fci + m2)i+a-i)= r(i + (/-i)£) 

The counterterm is then given, as usual, by minus the pole terms after expanding in e. After 
summing over the contributions of all webs at a given order and their counterterms for sub diver- 

2 

gences, all nonlocal terms (~ In -^) cancel as well as UV vertex counterterms and IR divergences, 
and we obtain the n-loop counterterm contributing at a; — > 1, which can be written as a series in 

" 1 (<^s{n')Win) 7dk' 



m=l 



„..y^, ,52) 







with purely numerical coefficients a^^. Because webs exponentiate, the counterterm for UV 
divergences in the perturbative expansion of a non-singlet parton distribution is given by 

oo m 1 

exp^EEltrl^l <^ 1^} (53) 

m=l n=l 
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in the limit x ^ 1, as indicated by the superscript A. Now it is trivial to extract the contribution 
to Pff. From ( ^41) and (|53|) we get 

Af^ = -na^l (54) 

As emphasized above, internal UV divergences, including the usual QCD divergences and 
divergences at the eikonal vertex, have to be renormalized. Further complications arise because 
coUinear/IR divergences cancel only after summing over all diagrams at a given order, so an 
individual diagram has in general UV singularities multiplying IR/coUinear singularities. Our 
method to resolve these technical problems is most transparent in light-cone ordered perturbation 
theory (LCOPT) |^T], 0, which is equivalent to performing all minus integrals of all loops. The 
Feynman rules for LCOPT are given in Appendix ^ for completeness. 

4.2 Summary of the Method in LCOPT 

Our method can be summarized as follows, details will be given below: 

1. We start with the expressions in LCOPT for the set of webs at a given order with a 
fixed coupling. The number of web-diagrams is much less than the number of all possible 
diagrams at a given order. Moreover, since we work in Feynman gauge, the number of 
possible webs is further reduced. For example, at order 2, as we will see below, only three 
diagrams contribute, aside from gluon self-energies. 

2. Ultraviolet divergent internal /c^ ^-integrals are regularized via dimensional regularization, 
with e > 0. At this stage we do not yet encounter IR/coUinear singularities since all 
integrals over transverse momenta are performed at fixed plus momenta. 

3. We add the necessary QCD counterterms and the counterterms for the eikonal vertex which 
has to be renormalized as a composite operator. As we showed in Section \i.2[ the sum of the 
latter cancels because of the recursive definition of webs and a Ward identity. However, in 
the intermediate stages the vertex counterterms are necessary to make individual diagrams 
UV finite. 

4. After elimination of the UV singularities we dimensionally continue to e < to regulate 
the IR/coUinear plus-integrals. It follows from the rules for LCOPT that all internal plus 
momenta are bounded by the total flowing into the minus eikonal. Therefore, the 
integrals over these internal plus momenta give no UV sub divergences. 

5. When we sum over the set of diagrams at a given order the IR divergent parts cancel, as 
well as the UV counterterms for the vertex, thus also the internally UV divergent vertex 
parts cancel. 

6. The final scaleless k± integral provides the UV counterterm contributing to the anomalous 
dimension (see Eqs. (Pl|)-([^) ). 
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The rules for LCOPT can be found in Appendix ^ I^T], they can be derived by performing 
the minus integrals. LCOPT is similar to old-fashioned, or time-ordered, perturbation theory, 
but ordered along the light-cone, rather than in In a LCOPT diagram all internal lines are 
on-shell, in contrast to a covariant Feynman diagram, which allows us to identify UV divergent 
loops in eikonal diagrams more easily. A covariant Feynman diagram is comprised of one or more 
LCOPT diagrams. 

We start by writing down all light-cone ordered diagrams of a given covariant Feynman dia- 
gram. All momenta in crossed gluon ladders have to be chosen independent of each other, such 
that they all flow through the eikonal vertex, since we seek the anomalous dimension for this 
vertex. Because ^ has no minus-component (cf. Eq. (^9])), we have g~ = in Eq. ([79| ) of the 



appendix when applying the Feynman rules for LCOPT in our case. This can be depicted graph- 
ically by contracting all propagators on the minus-eikonal (here jS) to a point, which coincides 
with the eikonal vertex. Two-loop examples can be found in Fig. |14[ Sometimes, numerators 
stemming from triple-gluon vertices or quark propagators cancel the corresponding propagators 
on the plus-eikonal ■ k = k~). Graphically, this can again be described by contracting these 
propagators to a point. Then we can read off easily from the various light-cone ordered diagrams 
the analytical expressions, whose k± i-mtegTals we perform, and renormalize. 

We need QCD counterterms and counterterms for the effective vertex. More specifically, by 
QCD counterterms we mean the usual gluon self-energy counterterms, as well as the counterterms 
for triple-gluon vertices and eikonal-gluon-eikonal vertices. The latter are UV divergent in any 
covariant diagram, however, this is not necessarily the case for all LCOPT diagrams found from 
a covariant diagram. Examples will be given below. Self-energies of the light-like eikonal lines 
vanish in Feynman gauge. Both types of counterterms are found via the (recursive) BPHZ- 
formalism, and the sub divergences are identified with the help of naive power-counting on a 
graph-by-graph basis. 

We will now illustrate our method by the rederivation of the 1- and 2-loop A-coefficients. 
4.3 Calculation of the 1- and 2-loop Coefficients Af\ Af 



The well-known |TT]] coefficients of the collinear parts of the splitting functions to one and two- 
loop order are given by 



A« = a, (55) 

(56) 



A^a'^ — 2^aK = -Ca 



I 18 6 9 



where Cq = Cp, Cg = Ca, Nf is the number of fermions, and Tr determines the normalization 
of the generators of fermion representation R, Tp = ^, Ta = Nc, with Nc the number of colors. 
We will now apply our method to the recalculation of these coefficients. 

The only web at order 1 in Feynman gauge is a single gluon exchanged between the two 
eikonal lines. The color weight is C*^^^ = Ca by definition (|37|) , and the web has no internal 
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momenta. Straightforwardly we obtain 



^i^wi^^c'-'^i^]!^] (i^rmr-^. (57) 




{2tt 

So at lowest order we get from Eq. (H) A'^J-^ = C^^^ = C^, as in (|55D . 



At order 2 we have the webs shown in Fig. |T3|, where we rotated the eikonal lines in the figure 
compared to the diagrams shown in Section |], to make the connection to Figs. |^ a) and ^ more 
evident. The original color factors are (compare to Fig. ^ 

cm = ^cp-^^Ca, 

CiW,) = C{W,) = -^Ca (58) 

for eikonal lines in the a-representation. The respective color weights of the webs c) and d) are 
the same as the original color factors, since they do not have decompositions (these diagrams are 
"maximally nonabelian"). The decomposition of diagram b) was shown as an example in Section 
^.1.2| , which resulted in a color weight 

CiW,) = -^Ca. (59) 

The contribution of Fig. [l^ a) is easily found from Eq. ( pTf ) and the well-known finite terms 
(see e.g. |20[) after renormalization of the of the gluon-self energy in the MS scheme, which is 
an example of what we called a QCD renormalization in the previous subsection: 

4'^''^^ = f^CACa + Y^C^Ca - iTnNfC, = - ^T«iV^) C,. (60) 



The first term in the first equality in Eq. (pOD stems from the gluon loop, the second term 
from the ghost loop. The last term is obviously the fermion loop contribution found from the 
expression for the graph. 



-w. 



(27r)" aa,a)Nf 



and its counterterm. 

The LCOPT diagrams obtained from the webs |13| b)-d) are shown in Fig. We see that 



due to the numerator {2k''' — k^) in the triple-gluon vertex, web d) contains two orderings on 
the light-cone; the factors of 2 and (—1) next to the eikonal vertices in the figure come from this 
numerator. Furthermore, for web b) it is important to route the momenta in the crossed ladder 
independently of each other, such that both of them flow through the vertex, to separate the 
subdivergence associated with the upper loop {k') from the overall UV divergence. 
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a) 



/3 



+ 






Figure 13: Webs contributing to (compare to Fig. |TT]): a) web of order 1 with 1-loop gluon 
self-energy inserted, b) the "crossed ladder" , c) and d) graphs with a triple gluon vertex. 



Now we determine the divergent 1-loop subgraphs for each web by naively counting the 
powers of transverse momentum components in numerators and denominators. The UV divergent 



subgraphs are marked with boxes in Fig. |14[ We see that for web c) and the first term of web d) 
we need a QCD counterterm for the triple-gluon vertex, whereas for the webs a) and d)(ii) we 
require vertex counterterms, as shown in the second column of Fig. |14[ Web d) (ii) is an example 
for a LCOPT graph with a triple-gluon vertex which does not need QCD renormalization, in 
contrast to loop-corrections to 3-gluon- vert ices in every covariant diagram. Due to the factor of 
(-1) in web d)(ii) the two vertex counterterms cancel each other, as announced above. The QCD 



counterterm, as shown in Fig. 15, is in the MS scheme for quark eikonal lines /3 given by 



c 



. 01 



1 /I 



IT 



(62) 



where g is the QCD coupling, = g^/i^in), and e > 0. 

The next step, after adding the appropriate counterterms to the respective graphs, is to 
perform the plus- momentum integrals. To do so, we dimensionally continue to e < 0, that is, to 
n > 4 dimensions. The results for the webs b)-d) and the counterterms for UV subdivergences, 
denoted by Z (omitting the vertex counterterms which cancel each other) are: 



-W, 



(2) 



(27r)""'^'^ 
dl'k 



-K J k+ 2 6 

/ \ / g 



(63) 



(27r)" 



W. 



(2) 



aa, c) 



dJ'k 



(2vr) 



-W, 



(2) 



aa, d) 



2, °fdk+ lT{2e) 



k+ A e 



X {B{1 - 6, -e) - 2B{1 -e,l- e)} 



(64) 
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Figure 14: LCOPT diagrams obtained from Fig. |T3|b)-d). The subgraphs in the dashed boxes are 
UV divergent 1-loop subgraphs, whose counterterms are shown in the second column. The grey 
boxes denote the eikonal vertex counterterms, whereas the grey blobs are the triple-gluon- vertex 



counterterms, shown in Fig. 15 
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2; (2^0^^^) = ^Jj2^nZ<im=C [-^j (4-) y ^2^(^Hi~^"47J• 

The color weight, as stated in Eqs. (|58D and (0), is C*'^'' = — for all diagrams. We notice 
that diagram d) gives the same contribution as its upside-down counterpart c), as expected, but 
only after adding different types of counterterms. 



a, /i a, /i 




Figure 15: QCD counterterm for the triple-gluon vertex, where PP denotes the pole part (omit- 
ting scheme-dependent constants). 



After summing over the contribution of the webs b)-d) and the counterterms, we see that the 
infrared poles l/{—e) in the Beta-functions cancel, as well as the vertex counterterms, leaving us 
with 

4^)'^)-'^) = ^a(2-^) (66) 

according to Eq. (Q). The contributions of all diagrams, ( |60D and (§^), result in the 2- loop 
coefficient (p6D, as announced. 



5 Higher Loops 

5.1 iV^-^-Terms in 

It is relatively straightforward to obtain a general formula for the A^j~^-contribution to the n- 
loop coefficient A^"^^ , since the only graphs involved are one- loop webs with n — m — 1 fermion 
bubbles and m counterterms for the fermion bubbles inserted into the gluon propagator. It is 
therefore a matter of simple combinatorics to obtain the a^NJ'^ contribution (compare to the 
one-loop expression Eq. (0) ): 



d"" _ onrrn-iMu^i dk^ 



(27r 
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n 



Ca{TRNfY~^-teim in A^^^ 



2 

27 

4 -ir + ^C(3) 

243 243 ^V"-"; 24|0 ^ 
6 ~729 ~ 729'^(^) ~ 4374 8T'^(^) 



Table 1: a,(/i2)"Ar"-i 



-contributions to the anomalous dimension Pjf. The expansion of 



is performed in terms of a^/Tr (cf. Eq. 



n-l 




n— m— 1 



2 \ (n-m)e 



m=0 



(67) 



r 1+ n 



[T{e)B{2-e,2-e)] 



n—m—l 



















47r /. 



The ^-pole in the expansion of the F- and Beta-functions in the sum is the contribution to 
the anomalous dimension (cf. Eq. (|5^) ). The contributions up to are given in Table |l|. 
They coincide with the corresponding values (the In iV-terms, or equivalently, the Si (iV)-terms) 
calculated by Gracey in fT^pl. 



5.2 Towards the Three-Loop Coefficient A 



(3) 
/ 



Vogt [0 obtained a numerical parametrization of the A^^'' from the known integer moments of 
the splitting function^: 



A 



(3) 



(13.81 ± 0.14) - (4.31 ± 0.02)TRNf 



27 



(68) 



We obtained the term proportional to Nj in the previous subsection, as listed in Table Now 
we will go on to compute the term proportional to Nf. All intermediate expressions are simple 
enough to be handled by the general algebraic computer program Mathematica |^^. For the 
calculation of the full A^^'' or even higher loops, however, an implementation of the algorithm 



into a more specialized computer algebra program such as FORM |34] may be desirable. 

The diagrams contributing to this term and their QCD counterterms are listed in Table 0, 
labelled in analogy to the two-loop case. We only have to compute the contributions from the 
Qf^i, part of the dressed gluon propagator, since the longitudinal parts 
Ward identity shown in Fig. |l^. This cancellation has been verified explicitly. 

The contributions to the set a) are easily computed to be 

1 



kfik^ cancel due to the 



A 



(3), a) 
/ 



-CATRNfCp. 



(69) 



^Note the different overall normalization of the anomalous dimension there. 

•^Note the expansion in (j^) there, whereas we expand in terms of (^) - see Eq. (||). 
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The contributions to the Nf-pait of the two-loop gluon self-energy inserted into a one-loop web 
(set g) ) give: 



/509 1 ^ / 55 



Ca{ + -C(3) ]-Cf{ + C(3) 

. V864 2^^ V V 48 ^ 



TRNfCp. (70) 



To compute the two-loop gluon self-energy, the occurring tensor integrals have been reduced to 
simple scalar one- and two-loop master integrals using the relations derived in We checked 
our calculations of the set g) against previous computations of the two-loop gluon self-energy 
in Feynman gauge, see for example |^6|. Note that this contribution has a term ~ Cp, which 



is not "maximally non-abelian" . The results of [36| include the longitudinal terms of the gluon 



propagator, which is dressed with a fermion bubble. Since these terms in the two-loop gluon 
self-energy, as stated above, cancel against the longitudinal parts in the remaining webs, Eq. 
(|70| ) does not contain these contributions. 

The expressions for the two-loop webs with a one-loop bubble-insertion are found easily from 
the corresponding two-loop expressions Eqs. (|63| ) -(|65| ) , taking into account the proper multiples 
of e in the Gamma- and Beta-functions due to the bubbles. The calculation of the triangles 
e) and f) is a bit more nontrivial. The resulting contributions can be found in the table. The 
results for e) and f) have been expanded in terms of e and Beta-functions using various identities 
tabulated in |[37[| . 

Since the infrared structure of the graphs is modified by the bubbles, which effectively raise 
the powers of the corresponding gluon propagators by e to a non-integer value, the upside-down 
counterparts do not give the same contributions. This asymmetry is not surprising, since we 
compute the coefficients collinear to the plus eikonal, thus introducing an asymmetry in how 
we treat the eikonal lines and the gluons attaching to them. However, we find that the sum of 
graphs in set d) gives the same contribution as the sum of graphs in set c), as can be seen from 
the tabulated expressions. 

The individual diagrams b)-f) have at most three UV (QCD) divergences and one IR/coUinear 
divergence, in addition to the overall scaleless ^"'"-integral. We observe that the diagrams with 
a one-loop counterterm for the fermion bubble and the one-loop counterterms for the triangle 
graphs have the same IR structure as the two-loop webs. Thus their IR divergences cancel 
separately from the rest of the diagrams. This implies that the collinear divergences have to 
cancel within the set of remaining diagrams, that is, within the set of webs with bubbles and 
the triangles. Moreover, we observe that the infrared divergences cancel within certain subsets 
of these graphs. Namely, they cancel separately between graphs b)(l), c)(l), and d)(l), between 
graphs b)(2), c)(2), and d)(2), as well as between c)(3), d)(3), e) and f). 

Summing over all contributions from graphs b)-f) we arrive at 

^ V288 54 3 ; ^ ^ 

We performed several checks of our computations. The infrared structure described above is 
one check of the results listed in Table Another check is the cancellation of non-local logarithms 
~ logM. Furthermore, the values of the and 1/e^-poles can be predicted from the one- 



33 



a^Nf has the following structure: 



and two- loop calculations performed in Section |4.3| [p8| . The sum of all diagrams contributing at 

(72) 



' * n. 



{2n 



n aa,Nf 



+ 



54 e^^ 
'167 

324 108 



Ca + -Cf 



Nf 



dk+ 
1+ 



The predictions of the higher poles in Eq. (|72D coincide with the poles obtained from the 
expansion of the calculated expressions listed in the table. 

Adding (|69D, ([70|), and (^) we obtain the term proportional to Nf contributing to the three- 
loop coefficient A^^^: 



A 



(3) 
Nf 



(209 Stt^ 7C(3)\ 

^^1216 



55 
48 



+ C(3) 



TrNjCp = -4.293 TnNjCp, (73) 



which agrees with the numerical prediction in Eq. (RHI 
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Table |: 



Web 


T^a ft or 




a) /7f^ 


2 


see Eq. (^) 


b)(l) 1/^ 

v 


2 


KM'^^^ft^B(2 £,2 + e) 


b)(2) ^ 


2 


irM3^^J^S(2 e,2 5)5(1 + 2£, 2e) 


b)(Cl) ^^^^ 


4 




0,1, ^ 


2 


+™34^j)(^|;+^2);5(2 .,2 .) 

X {5(1 - e, -e) - 25(1 - e, 1 - e)} 


c)(2) ^ 


2 


+irM3^^^f2")5(2 £,2 e) 

X {5(1 - e, -2£) - 25(1 - e, 1 - 2£)} 


0,3, ^ 


2 


+^3^^5(2 2 -e) 

X {5(1 - 2£, -e) - 25(1 - e, 1 - 2£)} 




2 


+fs:M3^^^(3^)5(2 £,2 e) 

X {^§5^5(1 - - 5(1 - 25, 1 - e)} 


d)(2, ^ 


2 


X {5(1 - £, -2e) - 4(ni±|li!5(i _ 1 _ 2£)} 


d)(3, ^ 


2 


+irM3^i^2^^|)5(2 £,2 e) 

X {5(1 - 2£, -e) - 25(1 - e, 1 - 2e)) 
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Continuation of Table 



Web 


Factor 


Contribution 


c)d)(Cl) ^ 


12 


_i^M2^^^iV, {5(1 - £, -£) - 25(1 - £, 1 - £)} 


c)d)(C2) Q 




4 


+fs:M2^''(2^')A/'£5(2 £,2 e) 








c)d)(C3) p 


?^ 


4 




c)d)(C4) 


8 


-'^^'^ 2 [l2 e 18 




2 


_irM3^^g^5(2 - e, 2 - e) {5(3 - 2£, -e) 
-5(1 - 2 - 26) + + (4C(3) - 2) e"] 




2 


KM^^^ft^B(2 e,2 e) {5(3 2^, e) 
-5(1 - 2 - 2£) - ^£ + (10C(3) - 2) 


e)f)(Cl) ^ 


4 


+i^M2-^miv^ {5(1 _ _e) - 25(1 - 1 - e)] 


e)f)(C2)^£J- 


4 
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Continuation of Table 0: 



Web 


Factor 


Contribution 




2 


see Eq. 



We introduced the following abbreviations: 







M 



(4vr) 



Ne = --In 



47r ' 



(74) 

(75) 
(76) 



Table 2: Webs contributing to the iVj-term of the three-loop coefficient A^^^ and their countert- 
erms (c.t.s), labelled (C). The cross denotes the counterterm for the fermion bubble. Similarly, 
the cross in the triple-gluon vertex denotes the counterterm for the fermion triangle. The grey 



blob represents the counterterm Fig. |T^, the grey blob with a cross is the 2-loop counterterm 
for the triple-gluon vertex with a fermion bubble inside. And finally, the black box denotes the 
fermion part of the 2-loop counterterm for the triple-gluon vertex. We omitted vertex countert- 
erms which cancel. We refrain from drawing all counterterms which give the same contribution. 
Instead, we indicate multiple contributions in the column "factor". A factor of 2 is due to the 
two complex conjugate contributions, and has already been taken into account in Eqs. (|69|) and 
fiTC 



37 



6 Conclusions 



We have developed and proved a method for the calculation of the coefficients proportional to 
of the non-singlet parton splitting functions, whose knowledge, for example, is important 



1- 

"or NNLL resummations. The method is based on the factorization properties of the splitting 
functions, and on the exponentiation of eikonal cross sections. 

We illustrated the method with the rederivation of the 1- and 2-loop coefficients A^^^ and 
A^'^\ as well as the NJ~^ terms at order n. We presented the result for the term proportional to 
Nf at three loops, which coincides with the approximate result obtained by Vogt [If. The full 
splitting functions at three loops are currently being computed by Moch, Vermaseren, and Vogt 
|T5[ with the help of the operator product expansion. Their results for the Nf-term provide 



a further check of our calculations. 

Although a calculation via the OPE provides the complete A^, or equivalently, x-dependence 
of the splitting functions, it involves a large number of diagrams and complex expressions at 
higher orders. A computation at three loops is a formidable task, and it seems unlikely that 
higher loop calculations will be completed in the near future. However, for certain observables, 
large logarithms due to soft and/or coUinear radiation become numerically important, and need 
to be resummed to as high a level in logarithms as possible. Our method, although limited to 
only the computation of A, has the advantage that higher-order computations are much less 
complex than within conventional methods, because the number of graphs is greatly reduced, 
and the expressions involved are relatively simple in LCOPT. Moreover, a fully computerized 
implementation of the algorithm should be straightforward. Therefore, the computation of the 
coefficients A at four or even higher loops may be within reach. 

A Rules for LCOPT 

As already stated above, these rules can be obtained by performing first all minus integrals 
1311,^. 



We start with forming all possible light-cone time orderings of a given covariant diagram. 

Only those configurations are kept which describe possible physical processes once the 
energy flow is specifled. 

For every loop we have a factor 

(27r)3-2= ^ ^ 

if we work in n = 4 — 25 dimensions. 
For every internal line we have a factor 



2V 



(78) 
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corresponding to the flow of plus momentum through the graph. 



• Every intermediate virtual state contributes a factor 



I 



(79) 



where we sum over all momenta comprising that virtual state, and where q is the external 
minus-momentum of the incoming state(s). 

• Every intermediate real state gives a momentum conserving delta-function: 



where the sum is over all momenta in that real state. 

• Since all lines are on-shell, we replace for every Fermion numerator its minus component 
by its on-shell value: 

r^^. (81) 
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